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I.  INTRODUCTION 


This  is  the  first  of  a  series  of  reports  to  be  prepared  in  support  of 
the  atmospheric  modeling  portion  of  a  comprehensive  Ground  Target  Signa¬ 
ture  (GTS)  modeling  program.  The  aim  of  the  atmospheric  modeling  pro¬ 
ject  is  to  develop  the  capability  of  predicting  atmospheric  effects  on 
ground  target  signatures  for  a  wide  range  of  meteorological  conditions. 
Accomplishment  of  this  objective  will  provide  the  information  required 
for  the  design  and  operation  of  terminal  homing,  target  acquisition, 
and  surveillance  systems.  It  will  also  influence  the  choice  of  operat¬ 
ing  wavelengths  for  artificial  sources  used  in  conjunction  with  these 
systems  CG.  Moreover,  it  will  provide  the  knowledqe  needed  to  discrim¬ 
inate  between  target  and  environmental  effects  EG. 

The  neodymium  (Nd5+)  laser  wavelength  of  1.06  pm  was  chosen  for  this 
initial  investigation  because  the  Nd^+  laser  is  currently  being  employed 
in  laser  guided  terminal  homing  missile  systems  [2],  The  purpose  of  this 
report  is  to  present  the  results  of  a  literature  search  (through  January 
1972)  and  theoretical  analysis  conducted  to  determine  the  characteris¬ 
tics  of  atmospheric  transmission  at  the  1.06  pm  wavelength.  The  last 
literature  survey  on  the  subject  was  conducted  by  Roy  and  Emmons  E3] 
in  1965,  but  analysis  of  atmosphere  effects  on  the  propagation  of  1.06 
pm  Nd3+  laser  radiation  has  been  treated  more  recently  E4-6U. 

Results  of  experimental  and  theoretical  transmission  studies  by  various 
investigators  are  used  to  estimate  the  transmission  characteristics  for 
1.06  pm  radiation.  The  atmospheric  effects  of  molecular  and  aerosol 
scattering  and  absorption  are  treated.  Other  effects,  such  as  turbulence, 
multiple  scattering,  nonlinearity,  and  quantum  mechanical  corrections 
are  neglected,  although  under  certain  experimental  conditions  these 
effects  may  be  significant  (e.g.,  transmission  through  a  heavy  cloud 
layer  where  multiple  scattering  becomes  significant  or  high  beam  inten¬ 
sity  where  nonlinear  effects  are  important).  Transmittance  models  which 
include  some  of  these  effects  at  various  wavelengths  including  1.06  pm 
will  be  treated  in  subsequent  reports  in  this  series. 

Section  2  is  a  discussion  of  Nd^+  laser  emission  and  a  quantitative 
description  of  atmospheric  effects  on  laser  radiation.  Problem  areas 
which  need  to  be  resolved  to  enable  one  to  make  reliable  transmission 
predictions  are  brought  out  in  this  section. 

Section  3  is  a  discussion  and  analysis  of  the  results  of  both  experimen¬ 
tal  and  theoretical  transmission  studies  made  by  various  investigators 
in  the  region  of  interest. 

Section  4  is  a  comparison  between  atmospheric  transmission  at  1.06  pm 
and  1.54  pm,  since  transmission  at  these  two  wavelengths  is  higher  than 
at  most  other  laser  wavelengths  available  [6]. 


Section  5  is  a  summary  of  the  report,  includes  some  concluding  remarks, 
and  provides  the  results  of  Section  3  in  tabular  form. 


2.  BACKGROUND 


2.1  Neodymium  Laser  Emission  Spectral  Characteristics 


Neodymium  (Nd)  lasers  are  receiving  much  attention  for  use  in  atmospheric 
laser  systems  because  of  their  emission  spectral  characteristics  that 
lie  in  a  region  for  which  the  atmosphere  is  relatively  transparent,  and 
because  they  can  be  operated  at  ^mbient  temperature  with  high  efficiency. 
The  trivalent  neodymium  ion  (Nd^+)  is  the  active  entity  of  these  lasers 
that  operate  at  or  near  1.06  pm  wavelength,  with  some  finite  linewidth. 

In  describing  their  spectral  characteristics,  it  is  necessary  to  differen¬ 
tiate  between  the  gross  and  fine  structure  of  the  emission  bands  [7], 

The  gross  structure  of  the  emission  bands  (the  generation  wavelength)  is 
dependent  on  the  host  material  and  on  the  concentration  of  the  doping 
material  [8].  The  fine  structure  of  the  emission  bands  is  mainly  depen¬ 
dent  on  the  resonator  properties  and  on  the  width  and  character  of  the 
luminescence  line  broadening  [7],  The  important  types  of  neodymium 
lasers  commonly  being  used  in  atmospheric  laser  systems  are  Nd:YAG  and 
Ndrglass.  These  two  lasers  are  briefly  discussed  below. 


An  Nd:YAG  laser  requires  relatively  low  threshold  pumping  power  and  can 
be  operated  in  either  a  continuous  wave  (CW)  or  pulsed  mode  with  reason¬ 
ably  high  powers.  Ytterlum  aluminum  garnet  (YAG)  as  a  host  material  is 
hard,  has  a  high  thermal  conductivity,  and  can  be  grown  with  good  optical 
quality  [8].  The  ^aser  generation  wavelength  is  at  1.065  pm  with  line- 
widths  of  about  I0A  [8].  Models  emitting  up  to  one  kilowatt  continuous 
power  at  1.06  pm  are  commerically  available  [9], 


Characteristics  of  Ndiglass  lasers  are  discussed  in  great  detail  in  the 
literature  [10-14],  The  major  disadvantage  of  glass  as  a  host  material 
is  its  low  thermal  conductivity  [10],  This  imposes  limitations  on  the 
diameter  that  can  be  used  for  CW  and  high  repetition  rate  operations  be¬ 
cause  the  glass  host  cannot  dissipate  heat  quickly.  An  important  factor 
which  strongly  influences  the  generation  spectrum  is  the  character  of 
the  luminescence  line  broadening  in  the  active  medium  [7],  For  an  in¬ 
homogeneous  I  y  broadened  luminescence,  the  width  of  the  laser  emission 
spectrum  can  be  quite  large.  Due  to  the  nature  of  the  glass  host,  line 
broadening  in  neodymium  glass  is  inhomogeneous  [7,10],  Snitzer  [15] 
reports  that  the  Nd^+  1.06  pm  glass  laser  emission  envelope  just  above 
threshold  is  approximately  3  x  I0”4  pm  wide  but  can  be  increased  to  approx¬ 
imately  1.3  x  I0"2  pm  by  increasing  the  pump  power.  He  further  reports 


that  individual 


10 


,-5 


laser  modes  of  widths  less  than 
anywhere  within  this  1.3  x  I0“2  pm  wide  emission  evelope. 


pm  can  be  obtained 
Hence,  the 


2 


spectral  composition  of  the  Nd:glass  laser  varies  within  a  single  flash 
[16],  Since  the  Nd:glass  laser  has  a  high  energy  storage  per  unit 
volume  and  the  glass  host  material  affords  considerable  flexibility  in 
size  and  shape,  very  high  output  powers  are  realizable  from  the  Ndtglass 
laser  [10,1 l], 

2.2  Transmittance 

In  this  report  the  simplified  model  [17]  of  an  atmosphere  composed  of  N 
plane  parallel  homogeneous  layers,  each  characteri  zed  by  an  extinction 
coefficient,  y.,  will  be  adopted  but  not  applied  to  the  1.06  pm  wave¬ 
length  becauseJof  lack  of  experimental  data  concerning  aerosol  and  mole¬ 
cular  extinction  coefficients  versus  altitude  at  this  wavelength.  The 
application  of  this  model  which  considerably  reduces  the  complexity  of 
the  problem  of  calculating  the  atmospheric  transmittance  of  radiation 
will  be  the  subject  of  subsequent  reports  in  this  series.  The  complexity 
of  the  problem  is  due  to  the  dependence  of  the  extinction  coefficient 
on  a  number  of  different  physical  properties  of  the  atmosphere  [18], 

In  essence,  this  plane  stratified  atmospheric  model  requires  that  the 
physical  properties  of  the  atmosphere  change  with  elevation,  but  not 
laterally  and  in  addition  requires  that  the  curvature  of  the  earth  be 
neglected. 

Consider  a  beam  of  monochromatic  collimated  laser  radiation  propagating 
along  the  path  A-£;  in  the  jth  layer  of  the  assumed  plane  stratified  at¬ 
mosphere.  The  transmittance  tj  is  given  by 

x  .  =  exp(-y .a£.)  ( I ) 

J  J  J 

where 

A  l.  =  AZ.sece.  (2) 

J  J 

Here  AZj  represents  the  increment  of  vertical  distance  and  0  is  the 
angle  between  the  vertical  and  the  direction  of  propagation.  The  total 
transmittance  of  this  model  atmosphere  is  given  by 


N 

t  =  n  t .  (3) 

j=i  J 
or 

N 

t  =  exp(-  Z  y ,a£  . ) .  (4) 

j=l  J  J 

The  extinction  coefficient  yj  can  be  expressed  as  the  sum  of  the  scat¬ 
tering  coefficient  aj  and  the  absorption  coefficient  kj,  i.e.. 


3 


(5) 


=  o .  + 
J 


The  quantities  c.-  and  kj  depend  on  the  kind  and  number  of  molecules  and 
aerosols  in  the  atmospheric  path  (see  section  2.6),  so  one  can  write 


and 


G  . 

J 


G  .  +  0  . 

"U  aj 


(6) 


where 


k. 

J 


k  ,  +  k  . 

mJ  aj 


(7) 


am  =  molecular  scattering  coefficient 

=  aerosol  scattering  coefficient 

k  =  molecular  absorption  coefficient 
m 

k  =  aerosol  absorption  coefficient. 
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All  of  the  above  coefficients,  y ,  a,  and  k,  are  defined  in  km  *  i  f  ht 
is  in  kilometers  (km). 


2.3  Molecular  Scattering 

The  molecular  scattering  (Rayleigh)  coefficient,  om,  as  a  function  of 
altitude  for  each  wavelength,  is  expressed  by 

a  (h,X)  =  aD(X)N(h)  (8) 

m  R 


where 


aR(X)  =  Rayleigh  cross  section  at  wavelength  X 
N(h)  =  molecular  number  density  at  altitude  h. 


Equation  (8)  is  used  to  compute  the  Rayleigh  coefficient,  am.,  character¬ 
izing  the  jth  layer  located  at  an  altitude  h.  The  Rayleigh  tiross  sections 
for  anisotropic  air  molecules  is  computed  from  Cl9] 


3  2  2 

m  -  8*1  (n  -  IP 

°R  X)  4  2 

3  X  N 


(9) 


where  f  =  3(2+6)/(6-76) ,  6  is  the  so-called  depolarization  factor,  n  is 
the  index  of  refraction  of  air,  and  N  is  the  molecular  number  density  at 
sea  level. 


4 


The  factor  f  is  included  in  Equation  (9)  to  compensate  for  the  degree 
of  molecular  anisotropy.  Gucker  and  Basu  [203  report  a  value  of  .035 
for  6. 


The  index  of  refraction  for  "standard  air,"  viz.,  dry  air  containing  0.03 
percent  by  volume  of  CO2  at  a  pressure  of  1013.25  mb  and  a  temperature 
of  288°K,  may  be  computed  from  the  dispersion  formula  of  Edlen  C2I3, 


(n-l)!06  =  64.328  +  2949?r  l.9-  +  255,40 

1 46- ( x”; 


4I-(X"2) 


(10) 


where  X  is  the  wavelength  in  micrometers.  In  terms  of  atmospheric  pressure 
P(mb)  and  temperature  T(°K)  the  refractive  modulus  (n-l)IO^  may  be  given 
by  a  less  precise  but  more  convenient  formula  C22],  i.e.. 


(n-l) I06 


77. 6P  , 

T  + 


0.584P 
T  X2 


(II) 


Equation  (II)  is  valid  for  wavelengths  from  0.2  pm  to  20  pm. 

2.4  Selective  Molecular  Absorption 

The  absorption  bands  of  the  atmospheric  gases  constitute  a  complex  super¬ 
position  of  many  spectral  lines;  therefore,  one  writes 

k  ,(v)  =  Ik  .  ,(v)  (12) 

mj  .  mi j 


where  km  j  j  is  the  molecular  absorption  coefficient  of  the  ith  spectral 
line  in  tne  jth  layer.  The  (i)  summation  is  over  all  molecular  absorption 
lines  of  all  absorbing  species  which  are  close  enough  to  the  wavenumber 
v  to  contribute  significantly  to  the  total  molecular  absorption  coefficient 


The  attenuation  of  laser  radiation  in  the  atmosphere  is  caused  in  part 
by  the  so-called  continuum  absorption  due  primarily  to  the  integrated 
effects  of  the  wings  of  nearby  strong  absorption  lines.  This  decomposi¬ 
tion  of  the  molecular  absorption  coefficient,  kmj(v),  may  be  expressed 
in  the  form 


k  . (  v)  =  Ek  _ .(  v)  + 


mj 


mrj 


Ik  .(v) 
c  mSJ 


(13) 


where  the  "r"  summation  is  over  all  spectral  lines  whose  centers  are 
within  the  spectral  width  Av  of  the  laser  radiation  and  the  "s"  summation 
extends  to  all  absorbing  species  beyond  the  limits  of  Av  which  affect 
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tne  transmission  in  Av.  The  first  term  is  the  total  selective  molecular 
absorption  coefficient.  The  second  term  is  the  continuum  molecular 
absorption  contribution  from  the  wings  of  the  many  distant  absorption 
lines  on  either  side  of  the  spectral  region  under  consideration.  The 
definition  of  the  term  "total  continuum  absorption  coefficient"  is  dis¬ 
cussed  in  the  following  section. 


The  molecular  absorption  coefficient  for  a  collision-broadened 

line  can  be  given  to  good  accuracy  in  the  lower  layers  of  the  atmosphere 
(altitude  <50  km)  for  the  pressure  range  from  about  lO^mb  to  about  10  mb 
by  [23] 


I 

- 3 - ? 

(v-v  .  .)  t  a,  . 
oij  ij 


(14) 


where 


S.  .  =  Ik  , .(v)dv. 

Ij  mij 


(15) 


Here  S|  ■  is  the  normalized  line  intensity,  a<  ,  is  the  half-width  at  half 
maximum,  v0jj  is  the  wavenumber  determining  tne  position  of  the  center 
of  the  ith  spectral  line  and  v  is  the  wavenumber  at  which  the  absorption 
coefficient  is  required.  The  Lorentz  line  shape  given  by  Equation  (14) 
is  valid  for  v  lying  within  a  few  tenths  of  a  cm- 1  of  the  line  center 
v0|j  [243.  It  has  been  shown  [243  that  when  (v-v0|)»aj  the  wings  of 
the  lines  decay  much  more  rapidly  than  Lorentz  lines.  Hence,  for  con¬ 
tinuum  molecular  absorption  calclatlons,  a  correction  factor  must  be 
introduced  into  the  Lorentz  shape  when  v  lies  beyond  a  few  tenths  of  a 
cm'*  from  the  line  center  vQj.  Line  shape  studies  are  needed  in  order 
to  find  the  correct  form  for  km(v)  that  would  be  valid  for  all  v. 


The  line  intensities  are  functions  of  temperature,  the  dependence  being 
different  for  different  spectral  lines.  If  S0  is  the  value  of  Intensity 
at  standard  temperature,  T0,  and  pressure,  PQ,  the  value  of  S  at  other 
temperatures  is  given  by  [253 


S(T,p)  =  Sop(TQ/T)aexp 


(16) 


where 


EM  =  energy  of  molecule  in  the  lower  state  of  the  transition 
responsible  for  the  spectral  line 
k  =  Boltzmann's  constant 

p  =  partial  pressure  of  absorbing  gas  (in  atmospheres) 

T  =  temperature  (in  degrees  Kelvin) 
a  =  constant  depending  on  absorbing  gas. 
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The  half-width  of  the  absorption  line  is  a  function  of  the  partial  pres¬ 
sure  or  the  absorbing  gas,  the  pressure  of  the  nonabsorbing  gas,  and  the 
temperature.  The  temperature  and  pressure  dependence  of  the  half-width 
is  approximately  [25] 


Here  aD  is  the  half-width  of  the  line  in  question  at  a  pressure  PQ  and 
a  temperature  Tc.  The  constant  b  depends  on  the  absorbing  gas.  The  ef¬ 
fective  pressure  Pe  is  usually  determined  from  the  empirical  relation 

[26  r 

P  =  P  +  ( B— I ) p  (  18) 

e 

where  P  is  the  total  atmospheric  pressure;  p  is  the  partial  pressure  of 
the  absorbing  gas;  and  8  is  the  self-broadening  factor  defined  as  the 
ratio  of  the  self-broadening  to  the  broadening  by  the  nonabsorbing  gas 
present.  For  small  amounts  of  absorbers  under  atmospheric  conditions 
P  >>  p  so  that  one  can  write  Pe%P. 

Calculations  of  the  integrated  transmittance  over  the  rotational  bands 
of  water  vapor  have  shown  that  the  temperature  dependence  of  this  quantity 
is  determined  mainly  by  the  temperature  dependence  of  the  line  intensities 
[27].  The  temperature  variation  of  the  line  intensity  may  change  by 
several  orders  of  magnitude  for  the  temperature  range  that  occurs  in  the 
atmosphere  [28].  At  1.06  pm  the  vibration- rotation  bands  dominate  absorp¬ 
tion  so  a  corresponding  study  of  the  effects  of  temperature  and  pressure 
on  the  half-widths  and  intensities  of  vibrational-rotational  lines  is 
needed  in  order  to  predict  accurately  what  these  effects  will  do  to  radia¬ 
tion  propagating  -through  a  real  atmosphere  under  any  weather  conditions. 

At  present  it  is  difficult  to  determine  experimentally  the  fine  structure 
of  the  vibrational-rotational  spectrum  (such  as  line  positions,  their 
intensities,  and  half-widths)  for  the  majority  of  lines  [29],  The  ex¬ 
perimental  difficulties  are  caused  by  the  abundance  of  absorption  lines, 
their  overlapping,  and  the  spectrometer’s  slit  function  which  makes  it 
impossible  to  determine  the  true  line  contours  and  masks  weak  narrow 
lines.  However,  line  strengths  and  half-widths  may,  in  principle,  be 
■  alculated  from  quantum  mechanics,  although  in  practice  the  calculations 
are  intractable.  An  attempt  to  predict  the  fine  structure  of  vibrational 
rotational  spectra  theoretically  has  been  made  by  Gates  et  al.  [30]  for 
the  2.7  band  of  H2O.  Although  the  authors  took  the  positions  of  the 
line  centers  from  experimental  data,  the  intensities  and  half-widths  were 
calculated  by  quantum  mechanical  methods.  However,  they  did  not  consider 
the  interaction  of  vibrational  and  rotational  motions  of  the  molecule. 
Failure  to  take  into  account  the  vibrational-rotation  interaction  in 
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water  vapor  may  result  in  errors  amounting  from  ten  to  several  hundred 
per  cent  C3lU.  Zuev  et  al.  C29,32j  have  taken  into  account  the  vibration- 
rotation  interaction  in  calculating  the  intensities,  half-widths,  and 
line  centers  by  quantum  mechanical  methods.  The  very  high  degree  of 
monochromaticity  of  laser  radiation  makes  this  degree  of  sophistication 
necessary  in  calculating  absorption  of  laser  radiation  in  the  atmosphere. 
This  is  true  for  the  neodymium  laser  whose  emission  evelope  can  be  rela¬ 
tively  broad  (as  wide  as  1.3  x  10“^  pm)  because  the  individual  laser  modes 
within  the  envelope  may  have  widths  less  than  10"^  ym  [ I 5D.  Presently, 
to  this  author's  knowledge,  no  quantum  mechanical  calculations  for  the 
molecular  absorption  coefficient  at  \  =  1.06  ym  for  any  absorbing  gas 
have  been  reported  in  the  literature. 

2.5  Total  Continuum  Absorption 

Practically  all  attempts  at  experimental  investigation  of  the  total  con¬ 
tinuum  absorption  have  been  made  in  the  8-12  ym  window.  Even  here,  the 
results  of  the  total  continuum  effect  are  rather  contradictory.  The 
absorption  coefficients  calculated  by  various  authors  for  the  8-12  ym 
spectral  region  differ  by  50  to  100$  [33].  One  of  the  major  problems 
is  that  the  definition  of  the  total  continuum  absorption  coefficient  is 
vague.  That  is  due  to  the  fact  that  attenuation  in  the  so  called  "atmos¬ 
pheric  windows"  caused  by  the  continuum  molecular  absorption  due  to  wind 
effects  of  nearby  strong  absorbing  lines  is  just  one  of  the  components 
of  the  overall  continuum  attenuation  of  the  radiation.  There  is  also 
attenuation  of  radiation  due  to  molecular  scattering,  aerosol  absorption, 
aerosol  scattering,  and  in  the  case  of  highly  monochromatic  lasers,  there 
may  be  strong  selective  molecular  absorption  by  atmospheric  gases  such 
as  H2O,  O2,  CO2,  and  O3.  Hence,  in  many  cases  all  of  these  attenuation 
factors  are  lumped  into  one  and  the  total  continuum  absorption  coefficient 
k^  in  the  spectral  width  Av  is  defined  by  the  equation  [33] 

exp{-kS£}  =  -~f  exp{-Ey ,A£.}dv  (19) 

T  fiviv  ,'j  j 

where 


t  =  Ihl. 

j  J 

is  the  total  path  and  the  extinction  coefficient,  yj,  includes  all  the 
following  terms: 


y .  =  a  .  +  £k  .t£k  .  t  o  .  +  k  ,. 
J  p  mrj  s  msj  aj  aj 


(20) 


Many  authors  simply  assume  that  all  but  the  third  term  of  Equation  (20) 
are  negligible  or  that  the  attenuation  effect  in  the  atmospheric  window 
Av  is  caused  solely  by  the  molecular  absorption  by  wings  of  nearby  lines. 
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The  continuum 
equivalent  to 
(13)) 


absorption  coefficient  k^  is  defined  in  this 
the  continuum  molecular  coefficient  given  by 


report  as 
(See  Equation 


k<? 

J 


The  total  continuum  absorption  coefficient  in  the  spectral  interval 
for  this  case  is  calculated  from 


exp{-kS£}  =  /  expf-Zk^A^ .}dv.  (21) 

T  4v  4v  r  ,  j  j 

Q 

The  quantity  k:  is  evaluated  at  a  position,  far  removed  from  the  centers 
of  the  absorbing  lines,  where  the  line  variations  with  wavenumbers  are 
small  in  the  interval  Av.  Hence,  the  summation 

Zk^A*: 

j  J 

may  be  treated  to  a  first  approximation  as  constant  in  the  spectral  in¬ 
terval  Av  so  that 


kS  =  i  ZkCA-e.. 

T  j  j 


(22) 


Q 

In  order  to  measure  kj  as  defined  by  Equation  (22),  the  effects  due  to 
the  other  attenuation  components  (see  Equation  (20))  must  be  considered. 


The  problem  of  the  temperature  dependence  of  ab'  'rption  in  the  continuum 
is  not  established  [33].  The  dependence  of  kc(  on  temperature  and 
pressure  should  be  resolved  if  predictions  of  c>.  Inuum  effects  are  to 
be  made  for  different  weather  conditions. 

2.6  Aerosol  Extinction:  Attenuation  of  Radiation  at  1.06  vm  Wavelength 
by  Dust,  Haze,  Fog,  Cloud,  Rain,  and  Snow 


The  transmission  of  near- inf rared  radiation  through  the  atmosphere  is 
largely  dependent  on  scattering  and  absorption  by  the  most  variable  and 
least  investigated  atmospheric  component,  the  atmospheric  aerosol.  The 
aerosol  may  consist  of  dust  and  combustion  products,  salt  particles,  in¬ 
dustrial  pollutants,  minute  living  organisms  and,  most  important,  water 
droplets. 

Mie  theory  can  be  used  to  derive  the  aerosol  absorption  (ka)  and  scat¬ 
tering  (oa)  coefficients  for  spherical  particles  (or  some  other  simple 
forms),  provided  the  complex  index  of  refraction,  number  density  and  size 
distribution  are  specified.  One  should  keep  in  mind,  however,  that  i rregu- 
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larities  in  particle  shape  will  disrupt  tne  resonance  inside  the  particle 
causing  marked  changes  in  the  scattering  [34].  Very  few  studies  have 
been  made  under  controlled  conditions  to  compare  the  calculated  and  ex¬ 
perimentally  measured  coefficients  of  radiation  attenuation  by  the  aerosol 
[35,36], 

An  assumption  which  is  often  made  is  that  absorption  by  the  scattering 
aerosol  particles  is  negligibly  small.  Since  over  most  of  the  spectrum 
aerosol  scattering  encountered  in  the  atmosphere  Is  dominant  over  aerosol 
absorption,  this  assumption  is  reasonable.  However,  there  are  exceptions 
where  the  scattering  particles  absorb  strongly  [37]  so  that  one  has  to 
be  careful  in  practice  before  making  this  assumption.  In  what  follows, 
the  aerosol  absorption  coefficient  will  be  assumed  negligible  and  thus 
the  terms  "aerosol  extinction  coefficient"  and  "aerosol  scattering  coef¬ 
ficient"  may  be  used  interchangeably.  Areas  where  the  magnitude  of 
aerosol  absorption  may  be  significant  will  be  indicated. 

The  atmospheric  aerosol  is  present  in  any  stratum  in  the  atmosphere  but 
with  a  highly  variable  concentration  both  in  space  and  time.  The  aerosol 
size  distribution  spans  a  large  range  of  particle  sizes  and  may  have  con¬ 
siderable  structure  as  found  by  Fenn  [38],  Goetz  [39],  Kondrat'yev  et  al. 
[40],  and  Rozenburg  [41].  Moreover,  many  of  the  natural  aerosol  particles 
which  are  mixtures  of  water-soluble  and  water-insoluble  components  will 
undergo  severe  size  fluctuations  as  the  relative  humidity  varies  [42,43], 
The  natural  aerosol  particles  vary  in  size  from  about  4  x  10“^  pm  radius 
to  about  10^  pm  radius  and  are  generally  classified  as  follows: 

a.  Aitken  particles  of  radius  n,  io  3  to  10  *  pm 

b.  Large  particles  of  radius  'v  10  *  to  1.0  pm 

c.  Giant  particles  of  radius  ^  1.0  pm. 

The  Aitken  particles  have  the  distinction  of  having  a  strongly  wavelength- 
dependent  absorption  spectrum  [41]. 

If  the  atmospheric  aerosol  is  treated  as  spherical  particles,  the  attenua¬ 
tion  of  radiation  by  aerosol  will  depend  on  the  particle  composition,  the 
particle  size  distribution,  and  the  number  density.  Consider  an  atmos¬ 
phere  where  the  aerosol  sizes  within  a  unit  volume  are  characterized  by 
a  size  distribution  function  f(r).  Various  size  distribution  functions 
have  been  offered  in  the  literature  to  describe  the  real  atmosphere  (e.g., 
Junge  [44],  Fenn  [38],  Fortzik  [45]).  The  aerosol  attenuation  coefficient 
yg  is  given  by 

N 

r2 

Y  (h,  n,  r.  A)  =  /  C_(n,  r,A)dN(r)  (23) 

a  b 
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where 


dN (  r)  =  N(h)f (  r)dr;  /  f(r)dr  =  I 


(24) 


and  dN(r) 


n 

A 

Cs 

rl'r2 


N(h) 


f  ( r) 

N  &  N 
ri|  r2 


number  of  particles  with  radii  between  r  and  r+dr  per  unit 
vo I ume 

aerosol  extinction  coefficient 

complex  index  of  refraction 
wavelength  of  incident  radiation 
cross  section  for  each  particle 

minimum  and  maximum  radii  of  the  aerosol  particles 

total  number  of  particles  with  radii  between  rj  and  per 
unit  volume  at  height  h 
size  distribution  function 

aerosol  number  density  limits  corresponding  to  the  limits 
r(  and  r2« 


So 


r2 

y  (h,n,r,A)  =  N(h)  /  Cc (n, r, A) f (r) dr. 

a  h 


(25) 


This  is  usually  written  in  terms  of  the  efficiency  factor  Qext  defined 
as  the  ratio  of  the  aerosol  cross  section  C5  to  the  geometrical  cross 
section  of  a  single  particle  nr?: 

0ex+  =  Cs(n,r,A)/*r  .  (26) 


Hence,  Equation  (24)  can  be  given  as 


Ya  =  N  /  -nr  Qex+(n,r,  A)f  (  r)dr.  , 


(27) 


From  Equation  (5/  it  is  seen  that  y  is  given  by 
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y  =  a  +  k 
a  a  a 


(28) 


where  each  of  the  coefficients  ya,  og,  and  ka  can  be  written  in  terms 
of  the  sum  of  the  corresponding  ith  particle  coefficients: 


(29) 


N 

=  Z  < 
i  =  l 


k 


a 


N 

Z  k 


i  =  l 


a 


These  definitions  of  ya,  aa,  and  ka  are  rigorous  only  if  multiple  scat¬ 
tering  is  ignored  C46l].  General  expressions  for  0  +  follow  from  Mie 
theory.  ex 

The  widely  varying  aerosol  characteristics  which  represent  a  specific 
atmospheric  condition  and  are  needed  for  the  theoretical  calculations 
of  the  extinction  coefficient  are  difficult  to  measure  in  practice. 
Hence,  the  most  effective  method  for  obtaining  the  spectral  variation 
of  the  aerosol  extinction  coefficients  is  by  direct  measurement  of  aero¬ 
sol  attenuation  of  radiation  at  various  wavelengths  and  for  different 
atmospheric  conditions. L47,48D. 


3.  PUBLISHED  DATA 


In  this  section  the  published  results  of  transmission  measurements  and 
theoretical  calculations  are  used  to  estimate  the  atmospheric  effects 
on  the  propagation  of  1.06  urn  laser  radiation.  The  only  available 
transmission  data  at  this  wavelength  was  obtained  with  apparatus  having 
limited  spectral  resolution  and  therefore  the  results  of  this  section 
are  first-order  approximations.  The  atmospheric  attenuating  mechanisms 
for  which  data  are  presented  are  molecular  scattering,  molecular  line 
absorption,  molecular  continuum  absorption,  and  attenuation  by  the  aero¬ 
sol.  Quantitative  estimates  of  the  effects  of  haze,  fog,  cloud,  and  rain 
are  provided.  Finally,  a  qualitative  estimate  of  attenuation  by  snowfall 
is  included. 


3.1  Molecular  Scattering  (o  ) 

Elterman  [I9[]  has  tabulated  the  theoretical  values  of  molecular  scatter¬ 
ing  coefficients  for  a  clear  standard  atmospheric  model  for  altitudes  up 
to  50  km.  Twenty-two  wavelengths  between  0.27  and  4.0  pm  are  considered. 
The  meteorological  range  at  sea  level  corresponds  to  about  25  km  at  0.55 
urn  wavelength.  The  meteorological  range  or  visual  range,  V,  is  that  dis¬ 
tance  for  which  the  transmittance  falls  to  2 %,  i.e.,  V  =  3.1912/g  where 
8  is  the  total  scattering  extinction  coefficient.  Figures  I  and  2  are 
plots  of  om  versus  \  and  ya  versus  X  based  upon  Elterman’s  data  [19], 

From  these  data  it  is  seen  that  the  molecular  scattering  coefficient  for 
1.06  pm  wavelength  radiation  should  be  of  the. order  of  8.2  x  lO-4  km" 
at  sea  level  and  approximately  6.1  x  I0"4  km  at  an  altitude  of  3  km. 


12 


SCATTERING  COEFFICIENT  (kn 


•awu'iMwww 


ATMOSPHERIC  SCATTERING  COEFFICIENTS  VERSUS 
WAVELENGTH  FOR  SEA  LEVEL  ALTITUDE 


WAVELENGTH  (/uni) 

V  Fiqure  I.  Rayleigh  and  aerosol  scattering  coefficients  at  sea  level  as  a 

function  of  wavelength  (from  data  by  Elterman  [I9J). 
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Figure  2.  Rayleigh  and  aerosol  scattering  coefficients  at  3  km  altitude  as  a 
function  of  wavelength  (from  data  by  Elterman  [19]). 
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These  are  of  the  same  order  as  the  values  obtained  for  the  molecular 
scattering  coefficient  crm  by  McClatchey  et  al.  [l?D  for  1.06  ym  radia¬ 
tion  along  a  path  for  five  different  atmospheric  models  (see  Table  I). 
it  will  be  seen  that  the  contribution  to  the  overall  attenuation  of  1.06 
ym  radiation  from  the  Rayleigh  component  is  small  and  may  be  ignored 
for  path  lengths  less  than  10^  km. 

3.2  Selective  Molecular  Absorption  (km) 

The  atmosphere  is  exceptionally  free  of  strong  absorption  lines  in  the 
1.06  ym  region  as  is  evident  from  solar  spectra  C49,50U.  Only  a  few  weak 
O2  and  H2O  lines  and  the  short-lived  LO2I2  and  DO2-N2I]  complexes  are 
active  in  this  region.  These  lines  are  difficult  to  observe  in  the  lab¬ 
oratory  even  with  high  resolution  and  long  paths  (^2  km). 

Figures  3-9  are  facsimiles  of  plates  26  to  29  from  the  Mohler  Solar  Atlas 
C49].  A  scale  can  be  arbitrarily  set  to  indicate  relative  absorption. 

No  observations  were  made  through  cloudy  or  hazy  skies.  An  upper  bound 
on  the  absorption  coefficient  may  be  obtained  for  the  maximum  selective 
line  absorption  coefficient  by  use  of0these  data.  The  strongest  absorp¬ 
tion  line  shown  is  located  at  10585.  IA.  If  one  assumes  an  equivalent 
16  km  sea- 1  eve  I  path  and  a  minimum  transmittance  of  35%  on  the  arbitrary 
scale  based  on  Taylor  and  Yates  data  [5lj,  Equation  (I)  yields  for  this 
line  (assuming  all  attenuation  is  due  to  line  absorption): 

k  ^  6.5  x  10  km  . 
m 

From  a  superficial  examination  of  the  spectrum,  it  appears  that  the 
average  baseline  transmittance  over  a  I0A  bandwidth  is  of  the  order  of 
95%.  This  yields  a  molecular  absorption  coefficient  for  this  region  of 
the  order  of 

.  base  ,  . «-3  .  - 1 

k  %  3  x  10  km  . 
m 

It  seems  reasonable  to  expect  the  line  absorption  contribution  to  the 
total  extinction  coefficient  to  fall  between  these  two  limits.  A  sta¬ 
tistical  analysis  by  Zuev  et  al.  C52]  of  atmospheric  transparency  to 
laser  radiation  at  various  wavelengths  reveals  that  the  absorption 
coefficient  k^  is  small,  not  exceeding 

k  ^  4  x  10  ^  km  * 

at  1.06  ym  under  the  conditons  of  Temperature  =  283°K  and  precipitable- 
centimeters  (pr-cm)  of  F^O  =  0.3  to  0.7. 
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TABLE  I 

Molecular  scattering  coefficients  reported  by  McClatchey  et  al.  [1 1 7H 
for  1.06  um  wavelength  for  five  different  atmospheric  models. 


Altitude 

(km) 

T  rop i ca 1 

a  ( km“  1 ) 
m 

Mid  latitude 

Summer 

a  ( km" 1 ) 
m 

Mi d latitude 

Winter 

o  ( km" ' ) 
m 

Subarctic 

Summer 

o  ( km“ ' ) 
m 

Subarctic 

Winter 

a  ( km" * ) 
m 

0 

8.04E-04 

8.20E-04 

8.9  IE-04 

8.38E-04 

9.39E-04 

0-1 

7.68E-04 

7.8IE-04 

8.43E-04 

7.98E-04 

8.77E-04 

1-2 

6.99E-04 

7.06E-04 

7.52E-04 

7.2IE-04 

7.70E-04 

2-3 

6.33E-04 

6.38E-04 

6.70E-04 

6.50E-04 

6.82E-04 
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Figure  5.  Atmospheric  absorption  lines  between  10,51 3. oX  and  I0,532.2A  (from  reference  C49D) 
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Figure  6.  A1w*pberlc  absorption  lines  between  I0.532.2A  ano  10.585.  |8  (from  reference  [49^3 
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Atmospheric  absorption  lines  between 


The  contribution  by  the  ZO2I2  and  CO2-N2D  comPlexes  has  been  under  in¬ 
vestigation  by  Di  anov-Klokov  et  a  I  .•  C  5  3l] .  They  have  shown  that  the  role 
of  the  band  of  the  complex  at  1.064  'Mm  is  dominant  over  the  role 

of  the  individual  O2  molecules.  Additional  data  on  the  bands  of  the 
complexes  in  the  near-infrared  are  needed  to  account  for  radiation  loss 
due  to  these  bands. 


3.3  Total  Cont  i  nuum  Absorpti on  OCj.) 


There  is  practically  no  experimental  data  concerning  the  water  vapor 
continuum  in  the  region  for  X<8ym  except  at  the  two  points  X  =  3.70  ym 
and  X  =  3.58  ym  where  k^  =  0.062  gm"  *  cm^  and  k*f  =  0.07  grrf  *  cm^, 
respectively  [543.  A  quantitative  assessment  by  Andreyev  and  Gal'tsev 
[33]  of  the  continuum  absorption  coefficient  for  the  region  X<8ym,  based 
on  the  numerical  analysis  of  the  fine  structure  data  of  the  water  vapor 
spectrum  presently  available  results  in 


kl.06  '  2  * 


10  "  gm  cm 


or 


or 


kl  .06  6 


10  molecules 


cm 


k^  06  ^  1.6  x  10  6  (atm-cm)  ^Tp. 

To  evaluate  the  contribution  of  k^  to  the  total  extinction  coefficient 
computed  from  Taylor  and  Yates'  data  [513  in  a  later  section  of  this  re¬ 
port,  consider  the  path  length  to  be  16  km,  the  temperature  to  be  278°K 
and  the  relative  humidity  (RH)  to  be  40$,  so  that  one  has 

c  -4-1 

kjv6x  10  km  . 


This  is  of  the  order  of  the  Rayleigh  coefficient  om,  although  at  305°K 
and  70$  relative  humidity,  k^  will  be  an  order  of  magnitude  greater. 

The  coefficient  of  continuum  absorption  by  H2O  for  Xn.  1.06  ym  computed 
by  Kozyrev  and  Bazhenov  [553  from  experimental  data  of  Elder  and  Strong 
[563  and  Streete  [573  on  the  attenuation  of  radiation  in  the  water  vapor 
windows  results  in 

,c  n  c  in-2  -I  2 

k|  06  ^  ^  ®  cm 


or 


1 .06 


n*  1.5  x 
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10  molecules 


cm 
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or 


c  —  5  “I 

k|  06  ^  4  x  *0  (a+m-cm)gTp. 

At  278°K  and  48$  relative  humidity  one  has 

k*?  'v  1 .5  x  10  2  km  * . 

I  .Ub 

This  value  extrapolated  from  available  experimental  data  is  25  times 
greater  than  the  theoretical  value. 

In  the  above  cases  it  is  not  clear  how  the  total  continuum  absorption 
coefficient  kj  is  defined.  This  may  be  the  reason  for  such  a  large  dis¬ 
crepancy  between  these  two  values.  The  largest  value  k^  =  1.5  x  I0“^km-* 
seems  somewhat  of  an  overestimate  but  is  useful  as  an  upper  limit  for  k^ 
at  T  =  278°K  and  RH  =  48$.  At  T  =  29I°K  and  RH  =  55$,  this  limit  increases 
to  k^  =  4  x  I0“^km”l.  Thus  it  is  seen  that  k^  is  a  function  of  temperature 
and  relative  humidity. 

3.4  Aerosol  Extinction:  Attenuation  of  Radiation  at  1.06  pm  Wavelength  by 
Dust.  Haze,  Foq,  Cloud,  Rain,  and  Snow 

McCormick,  Lawrence,  and  Crownfield  [58]  calculated  backscatter i ng  and 
total  cross  sections  for  a  Junge  [44]  size  distribution  f(r)  of  spherical 
aerosol  particles  with  a  refractive  index  of  1.5.  These  calculations  were 
made  for  incident  laser  radiation  at  wavelengths  of  0.3472  pm,  0.5300  pm, 
0.6943  pm,  and  1.06  pm,  for  different  values  of  the  size  distribution 
parameter  where 

f(r)  =  br  ^v+^.  b  =  constant 

and  for  four  sets  of  particle  radii  limits  ( r  | ,  r2>.  Only  the  aerosol 
component  was  considered  in  their  computations.  For  v  =  3  and  (i  j  =  0.04 
pm,  r2  =  10.0  jpm),  their  results  for  the  average  scattering  cross  sections 
per  particle,  Cs,  at  the  different  wavelengths  are 

C  (A  =  0.3472  pm)  =  2  x  I0_l0cm2 

— S  -  1 0  2 

C  (A  =  0.5300  pm)  =  1.32  x  10  cm 

— s  -  1 0  2 

Cs  (A  =  0.6943  pm)  =  1 .01  x  10  'crn 

Cs  (A  =  1.06  pm)  =  6.56  x  10  cm2  .  • 

If  a  "clear”  atmosphere  (23  km  visibility),  is  assumed,  the  particle  density 
at  sea  level  is  given  by  McClatchey  et  al.  Cl7]  approximately  as  2.8  x  I03 
particles  cm"^  so  that  the  average  scattering  coefficients  computed  from 
McCormick  et  al.  [58]  data  are 
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o  (A  =  0.3472  pm) 

_a 

a  (A  =  0.5300  pm) 
a 

a  (A  =  0.6943  pm) 
0 

a  ( A  =  i .06  pm) 

9 


"v  5,6  x  10  2km  * 
-v  3.7  x  IO-2km"* 1 
'v-  2.8  x  I0""2km  * 
-v  1.8  x  l(f2km"  1 


It  is  seen  that  the  average  scattering  coefficient  for  A  =  1.06  pm  is 
smaller  than  the  scattering  coefficients  at  the  shorter  wavelengths.  It 
may  be  additionally  inferred  that  if  aerosol  scattering  is  the  dominant 
mechanism  for  reducing  the  transmission,  the  1.54  pm  wavelength  should 
have  a  higher  transmission. 


Curcio,  Knestrick,  Cosden,  and  Drummeter  [59]  have  reported  on  the  atten¬ 
uation  of  radiation  by  aerosol  scattering.  They  reached  the  conclusion 
that  the  effective  particle-size  distribution  for  a  particular  day  at  any 
location  near  the  coast  can  generally  be  approximated  by  a  combination 
of  continental  and  maritime  distributions  [603.  A  typical  measured  rela¬ 
tionship  between  the  atmospheric  aerosol  extinction  coefficient  ya  and 
the  wavelength  for  a  particular  day  in  the  Chesapeake  Bay  area  is  shown 
in  their  paper  [593.  They  corrected  the  data  to  show  only  the  attenua¬ 
tion  due  to  The  atmospheric  aerosol.  From  their  data,  the  aerosol  extinc¬ 
tion  coefficient  at  A  =  1.06  pm  for  a  16  km  path  and  25  km  visibility  is 
approximately  given  by 


.  T  =  29I°K;  RH  =  55% 
y  s.  |.|3  x  I0"1  km-1 

pr-cm  of  H^O  =  4. 

From  Elterman's  data  [193,  it  is  seen  (Figures  I  and  2)  that  in  the  1.06 
pm  wavelength  region,  the  Rayleigh  scattering  coefficient  am  is  of  no 
importance  as  compared  with  the  aerosol  extinction  coefficient  ya.  The 
value  for  y_. ( A  -  1.06  pm)  as  qiyen  by  Elterman  (for  meteorological  range 
-  25  km)  is  va  ■>.  1.13  x  IO-fkm-', 

Following  is  a  discussion  of  results  concerning  the  attenuation  of  1.06 
pm  wavelength  radiation  by  important  types  of  aerosols:  humid  haze,  fog, 
cloud,  rain,  and  snow.  Experimental  measurements  of  attenuation  by  dry 
dust  and  other  aerosols  not  included  in  the  above  types  are  not  presented 
because  of  lack  of  data.  Realistic  theoretical  estimates  of  dry  particles 
such  as  these  (dust,  industrial  pollutants,  etc.)  cannot  be  made  because 
these  particles  are  very  irregular  in  shape.  Irregular  particles,  as 
mentioned  in  Section  2.6,  do  not  scatter  radiation  in  entirely  the  same 
way  that  spherical  particles  do  and  hence  Mie  theory  is  not  applicable. 

The  water  envelope  on  aerosols  tends  to  make  the  spherical  shape  assump¬ 
tion  more  realistic,  although  how  realistic  is  open  to  question. 
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Haze  is  defined  as  that  finely  dispersed  fraction  of  the  atmospheric 
aerosol  with  particle  radius  sizes  lying  between  0.1  pm  to  0.6  pm.  Hazes 
are  formed  when  the  relative  humidity  is  less  than  80$.  When  the  relative 
humidity  Increases  to  about  80$  or  greater,  a  haze  Is  replaced  by  a  foggy 
haze,  which  in  turn  is  replaced  by  a  fog  as  the  saturation  point  is  reached. 
The  essential  characteristic  of  a  fog  as  well  as  for  a  foggy  haze  is  the 
existence,  against  the  background  of  a  continuous  particle  size  distribution, 
of  two  stable  narrow  ranges  of  average  sized  particles  [48].  For  foggy 
hazes  the  ranges  of  particle  radii  are  I  to  5  pm  and  12  to  15  pm,  while 
for  fogs,  they  are  8  to  12  pm  and  18  to  25  pm.  The  radii  of  the  droplets 
that  dominate  the  scattering  characteristics  by  clouds  are  between  5  and 
20  pm.  The  average  radii  for  rain  droplets  exceed  100  pm. 

Approximately  90  to  95$  of  the  aerosol  formations  characterizing  the 
state  of  the  atmosphere  under  real  meteorological  conditions  are  hazes 
[48].  Rosenburg  £48}  reports  that  a  Junge  size  distribution  should  not 
be  used  to  calculate  extinction  coefficients  in  the  case  of  humid  hazes, 
because  of  the  water  envelope.  According  to  Gebble  et  al.  H6 1 H  and  Sin¬ 
clair  [62],  natural  haze  is  approximately  transparent  in  visible  and  near- 
infrared  regions,  and  according  to  Arnulf  et  al.  [63]  the  transmission 
increases  markedly  with  increasing  wavelength,  from  the  visible  to  10  pm. 
Arnulf  et  al.  measured  the  transmittance  t  through  haze  by  use  of  the 
expression 


where  d  is  the  optical  density  per  unit  length  and  x  is  the  atmospheric 
path  length.  The  optical  densities  and  extinction  coefficients  obtained 
by  Arnulf  et  al.  for  the  case  of  haze  are  shown  in  Figure  10.  The  optical 
densities,  d,  are  related  to  the  extinction  coefficients  by 


Y  =  2.3d. 

At  1.06  pm  the  densities  fall  in  the  interval  I  km  *  >  d  >  0.5  km  *.  The 
extinction  coefficient,  Y^aze*  *n  '*’h's  case  **es  1°  +he  Interval 

2 . 3  km  *  >  Yu  >  i .  1 5  km  * 
haze 

indicating  that  for  a  1  kilometer  path,  the  transmittance  t  falls  in  the 
interval 


10$  <  x  <  32$. 

The  value  for  y  computed  by  Arnulf  et  al.  probably  also  includes  absorp¬ 
tion  by  the  haze,  because  it  has  been  shown  by  Robinson  [64]  and  Kondratyev 
[65]  that  it  is  likely  that  haze  not  only  scatters,  but  also  absorbs  sig¬ 
nificantly  in  the  visible  and  infrared  spectral  regions. 
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Figure  10.  Scattering  coefficients  and  optical  densities  for  the 
case  of  haze  (from  reference  C&33). 
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Arnulf  et  al.  [63]  also  investigated  the  extinction  of  radiation  by  fogs. 
Figure  II  gives  the  extinction  coefficients  and  optical  densities  of  sev¬ 
eral  types  of  fogs  as  compiled  by  them.  As  seen  from  the  figure,  a  partic¬ 
ular  stable  fog  exhibits  approximately  the  same  opacity  at  all  wavelengths 
shown  out  to  4  ym  at  which  point  it  becomes  more  transparent.  Moreover, 
the  transmittance  of  the  various  fogs  at  10  um  Is  observed  to  be  higher 
than  at  all  other  wavelengths  shown  In  the  figure.  It  should  be  noted 
that  the  mean  radii  of  the  droplets  of  most  fog  Investigated  by  Arnulf 
et  al.  fall  in  the  range  2.3  to  3.0  ym,  which  is  substantially  less  than 
the  results  of  other  investigators.  The  maximum  measured  drop  radius 
never  exceeded  15.  ym. 

Kurnick,  Zitter,  and  Williams  [66]  measured  the  transmittance  of  natural 
fogs  at  various  wavelengths  and  found  a  rronotonic  decrease  of  the  extinc¬ 
tion  coefficient  with  wavelength. 

The  controversy  as  to  what  degree  the  optical  properties  of  fogs  determined 
experimentally  with  special  chambers  containing  artificial  fogs  can  be  used 
to  predict  transmission  through  natural  fogs  has  not  yet  been  resolved 
[67,68]. 

Gates  and  Shaw  [69]  measured  transmission  through  various  cloud  types  using 
the  sun  as  a  source.  They  report  that  transmission  through  clouds  is 
always  slightly  higher  in  the  8.0  to  12.0  ym  wavelength  region  than  in 
the  0.48  to  5  ym  wavelength  region.  Since  calculations  of  pure  scattering 
by  clouds  show  better  transmission  at  the  shorter  wavelengths  than  at 
the  longer  wavelengths,  the  authors  conclude  that  the  water  vapor  bands 
in  the  near  infrared  absorbed  sufficiently  to  reverse  the  trend. 

Carrier,  Cato,  and  Von  Essen  [70]  computed  extinction  coefficients  for 
eight  specific  cloud  models  using  the  exact  Mie  theory  for  incident  ra¬ 
diation  at  wavelengths  of  0.488,  0.694,  1.06,  4.0  and  10.6  ym.  Their  re¬ 
sults  showed  that  there  is  no  clear  advantage  of  one  wavelength  over  another 
for  improving  through-cloud  transmission,  although  backscattering  is  mark¬ 
edly  reduced  at  the  longer  wavelengths.  The  largest  source  of  uncertainty 
of  their  results  is  the  neglect  of  multiple  scattering.  Table  II  is  a 
summary  of  calculated  extinction  coefficients  versus  wavelength  for  the 
eight  models  considered.  Wide  variations  in  these  values  can  be  expected 
for  cloud  models  that  differ  from  those  used  in  this  investigation. 

Scattering  by  rain  is  a  multiple  scattering  phenomenon  which  is  difficult 
to  analyze  quantitatively.  However,  an  estimate  of  attenuation  by  rain 
can  be  made  by  assuming  single  scattering  theory  and  by  taking  notice  that 
multiple  scattering  effects  will  only  lower  the  measured  extinction  coef¬ 
ficient.  Thus,  by  Equation  (27),  i.e., 

r2 

Vain  =  N  1  *r2Qex+(n,r,X)f(r)dr  (30) 

rl 
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Figure  II.  Scattering  coefficients  and  opt! 
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TABLE  I  I 


Extinction  coefficients  versus  wavelength  for  eight  specific  cloud  types 
(from  reference  [70]). 

EXTINCTION  COEFFICIENT  y  (km"1) 

Wavelength 


Cloud  Type 

0.488  urn 

0.694  um 

1.06  pm 

4.0  pm 

10.6  pm 

Nimbostratus 

128 

130 

132 

147 

136 

Altostratus 

108 

109 

1  12 

130 

83.9 

Stratus  1 1 

100 

101 

103 

114 

104 

Cumulus  Congestus 

69.2 

69.8 

71.3 

81 .0 

67.6 

Stratus  1 

66.9 

67.9 

69.7 

90.  1 

42.8 

Cumulonimbus 

43.5 

43.8 

44.4 

48.2 

50.9 

Stratocumu 1  us 

45.3 

46.0 

47.1 

59.6 

24.8 

Fai r  Wx  Cumulus 

21.0 

21.3 

21.9 

27.6 

11.7 
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an  approximate  upper  bound  for  the  extinction  coefficient  for  rain  can 
be  found  for  a  given  concentration  N.  Now  at  wavelengths  between  0.4  pm 
and  15  pm,  the  rain  drops  can  be  considered  as  large  droplets  with  Qext 
virtually  equal  to  two. 


Hence, 


Here 


y  ,  =  2N  /  nr2f(r)dr  =  2G. 

'rain 


G  =  N  /  irr2f  ( r)dr 


(31) 


(32) 


is  the  geometrical  cross  section  of  drops  In  unit  volume.  Thus  the  ex¬ 
tinction  coefficient  for  rain  given  by  Equation  (31)  is  independent  of 
wavelength  between  0.4  pm  and  15  pm. 


Polyakova  CV I U  has  shown  from  experimental  studies  that  the  magnitude  of 
Train  ^rom  0,4  to  15  pm  wavelength  can  be  obtained  from  rainfall  intensity 
data.  Polyakova  established  the  empirical  relationship 


Yrain=0*2IJ 


0.74 


(33) 


where  J  is  the  rainfall  rate  in  mm/hr  and  yrajn  is  the  magnitude  of  the 
extinction  coefficient  for  rain  in  units  of  km-1.  Figure  12  gives  rain¬ 
fall  rate  versus  attenuation  by  rain  for  near  infrared  wavelengths. 

The  extinction  coefficient  for  rain  can  thus  be  estimated  from  Equation 
(31)  if  the  analytical  form  of  the  drop  size  distribution  and  the  drop 
number  density  are  known  or,  from  Equation  (33),  if  the  rainfall  intensity 
is  known. 

The  value  for  Ypain  ^  =  1*06  pm)  calculated  by  Equation  (33)  for  I  mm/hr 
rate  of  rai nfal I  is 

Y ra j n ( X  =  1.06  pm)  =  2.1  x  10  *  km  * 

and  for  10  mm/hr  rate  of  rainfall 

Y  .  (X  =  I .06  pm)  =  1.15  km  * . 
rain 


Snow  consists  of  scatterers  of  very  complex  form;  this  makes  theoretical 
predictions  of  attenuation  at  different  wavelengths  extremely  difficult. 
Although  the  extinction  coefficient  for  rain  is  independent  of  wavelength 
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TOTAL  EXTINCTION  BY  RAIN  FOR  ALL  LASER 
WAVELENGTHS  BETWEEN  0.4  AND  150  Mm 


RAIN  RATE  (mm/hr) 


Figure  12,  Rainfall  rate  (mm/hr)  versus  attenuation  (dB/km). 
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between  0.4  urn  and  15  um,  the  same  is  not  true  of  snow.  Wilson  and  Penzias 
[72]  have  observed  the  wavelength  dependence  of  the  measured  attenuation 
coefficient  in  snowfall  where  they  found  that  the  attenuation  is  greater 
at  10  um  wavelength  than  at  0.63  um  or  3.5  um  wavelength.  Measurements 
by  Kabanov  and  Pkhalagov  [73]  indicate  that  attenuation  by  snowfall  at 
10  um  wavelength  is  greater  than  at  I  um  wavelength.  (Attenuation  is 
still  severe  at  I  um  wavelength,  about  between  that  due  to  fog  and  that 
due  to  rain  for  the  same  amount  of  equivalent  liquid  water  [74].)  This 
suggests  that  transmission  through  snowfall  is  better  at  1.06  um  (Nd^+laser) 
than  at  10.6  um  ((X>2  laser)  wavelength.  This  inference  has  not  yet  been 
experimentally  tested. 

There  are  few  experimental  data  available  on  the  aerosol  attenuation 
coefficient  ya(h)  in  the  infrared  at  various  altitudes  at  the  presenT 
time  [75]. 


4.  COMPARISON  BETWEEN  TRANSMISSION  AT  1.06  um  AND  1.54  um 

Berlinguette  and  Tate  [76]  have  made  spectroscopic  measurements  of  abso¬ 
lute  transmission  for  atmospheric  paths  for  the  wavelength  region  0.9  um 
to  5.9  um  with  horizontal  path  lengths  of  0.67  km  and  1.2  km.  Selective 
window  transmissions,  in  terms  of  precipltable  water  vapor,  were  established 
for  the  0.92  um  to  1 . 10  um  window  and  other  spectral  windows  for  path 
lengths  of  0.007  km,  0.01  km,  0.27  km,  0.67  km,  and  1.2  km.  The  amounts 
of  precipitable  water  in  the  paths  considered  varied  from  0.02  to  7.4  mm. 
Unfortunately,  errors  of  +3JK  appeared  to  be  involved  in  the  calibration 
of  their  apparatus  so  thaT  little  reliance  can  be  placed  on  the  absolute 
scale  of  their  transmission  plots.  One  can,  however,  compare  the  value 
of  the  extinction  coefficients  for  1.06  um  neodymium  (Nd^-)  and  1.54  um 
erbium  (Er^+)  laser  radiation  from  their  data.  For  the  longest  path  (1.2 
km)  with  resolving  power  of  about  200  at  1.06  um  wavelength  and  7.4  pre¬ 
cipitable  millimeters  of  water  vapor  in  the  path,  their  data  show  that 

y |  Qg  %  1 . 1 3  x  10  *  km  * 

-2  -I 

Yl  54  ^  4*3  x  10  km 

That  is, 

*1.06  '  2>6»l.54- 

In  terms  of  precipitable  mm  of  h^O: 

Y|  Qg  n,  1.8  x  10  2  mm  '(l-^O) 
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compared  to  the  erbium  case  (atmospheric  transmission  characteristics  of 
Er^+  discussed  in  reference  C77]) 

yI.54  'v-  6.9  x  I0"3  mm_,(H20). 

Taylor  and  Yates  C5I]  have  made  infrared  transmission  measurements  from 
0.5  pm  to  15  pm  wavelength  over  three  horizontal  sea- 1  eve  I  paths  of  0.305 
km,  5.5  km,  and  16.25  km  at  Chesapeake  Bay.  Figure  13  displays  their 
results  for  a  relatively  clear  winter  day  under  the  conditions  shown. 

They  covered  the  band  0.5  pm  to  1.8  pm  with  a  resolving  power  of  353  at 
1.06  pm.  The  extinction  coefficients  calculated  from  their  data  for  both 
Nd^+  and  Er^+  laser  wavelengths  are 

a)  5.5  km  path:  Y|  q6  'v  1 . 15  x  10  *  km  * 

Y |  ^  ^  7.8  x  10  2  km  * 

or  in  terms  of  precipitable  water  vapor  in  the  path* 


yI.06  *  4*63  X  ,0"2 

-1 

mm 

<h2o> 

-2 

-| 

(h2o> 

Yl.54  *  3J4  x  10 

mm 

-2 

b)  16.25  km  path: 

Yl  .06 

'v  6.3  x 

10  z 

km 

-2 

Yl  .54 

'v  4.7  x 

10 

km 

or  in  terms  of  precipitable  water  vapor  in  the  path: 


.06  '  2‘°  * 

,_-2  "I 

10  mm 

<h2o> 

-2  -1 

.54  '  L5  x 

10  mm 

<h2o> 

R.  K.  Long  [6]  has  examined  transmission  measurements  made  by  Streete  [57], 
At  medium  resolution,  the  data  indicated  that  transmission  at  1.06  pm  may 
be  higher  than  that  at  1.54  pm.  This  i s  at  odds  with  the  results  of  this 
report  which  are  based  essentially  on  Taylor  and  Yates  data.  However,  at 
higher  resolution,  Streete’s  data  indicated  that  the  absolute  transmittances 
are  comparable  and  transmission  at  1.536  and  1.544  pm  may  be  higher.  The 
conclusion  of  the  author  corroborates  that  of  Long  [6]  that  Streete’s 
data  are  not  any  more  accurate  than  those  of  Taylor  and  Yates.  Thus  the 
problem  regarding  the  transmission  at  1.0b  pm  vs  1.54  pm  is  still  unre¬ 
solved.  In  any  case,  it  is  safe  to  say,  transmission  at  these  two  wave¬ 
lengths  is  higher  than  at  most  of  the  other  laser  wavelengths  available. 

The  problem  of  which  wavelength  has  a  higher  transmission  awaits  higher 
resolution  experimental  studies  at  these  two  wavelengths. 


35 


8  g  s  s  a  ° 

U.N30d3d)  NOISSIWSNVdl 


w 

Q) 

+- 

<0 

>- 

"D 

C 

ID 

L. 

o 

> 

10 

I- 


(/)  I — I 
c 

10  <D 

i-  u 
+-  c 


—  (0 
L  4- 

Q)  Q) 
x:  L- 
a. 


<  ~ 


<D 

L. 

D 

O) 

Ll. 


36 


WAVELENGTH  (MICRONS) 


5.  SUMMARY  ANC  CONCLUDING  REMARKS 


The  scattering  mechanism  is  of  considerable  importance  in  the  propagation 
of  1.06  pm  wavelength  radiation  because  this  wavelength  lies  in  the  0.92 
to  1.10  pm  "atmospheric  window"  of  weak  molecular  absorption.  The  usual 
assumption  made  under  these  circumstances  is  that  the  total  extinction 
coefficient  measured  in  this  region  of  weak  molecular  absorption  is  the 
total  scattering  coefficient  (Rayleigh  plus  Mie).  The  contributions  by 
the  continuum  absorption,  aerosol  absorption,  and  selective  line  absorp¬ 
tion  are  thus  assumed  negligible.  This  atmospheric  window  nonetheless 
is  filled  with  numerous  weak  lines  [O2,  H2O,  [02^2#  [O2-N2]]  whose  central 
frequencies  readily  absorb  laser  energy,  the  aerosol  may  absorb  signifi¬ 
cantly  and  the  magnitude  of  the  cortinuum  is  still  undecided.  For  these 
reasons  the  assumption  that  the  total  extinction  coefficient  measured  in 
this  window  region  is  the  total  scattering  coefficient  is  unreasonable 
for  the  case  of  atmospheric  attenuation  of  1.06  pm  laser  radiation.  To 
resolve  this  problem  (l.e.,  account  for  the  various  attenuation  components 
in  the  window  regions)  the  characteristics  of  the  aerosol  must  be  speci¬ 
fied,  the  lines  in  the  atmospheric  spectra  must  be  specified  to  no  less 
rhan  0.1  cm“l  accuracy  and  the  strength  of  the  continuum  must  be  established. 
For  lower  accuracy  of  line  positions,  errors  in  the  determination  of  the 
absorption  coefficient  may  reach  tens  and  hundreds  of  percent  and  even 
higher  [31].  For  low  visibi I ities  (<  10  km),  however,  the  aerosol  scat¬ 
tering  mechanism  will  dominate  the  attenuation  of  1.06  pm  electromagnetic, 
radiation  and,  hence,  molecular  scattering  and  absorption  losses  will  be 
of  minor  significance. 

It  should  be  noted  that  selective  window  transmission  calculations  by 
Berlinguette  and  Tate  [76]  were  not  subjected  to  the  calibration  error. 

Their  data  showed  that  for  amounts  of  precipitable  water  less  than  7.4 
mm  in  the  path,  the  transmission  in  the  0.92  pm  to  1.10  pm  window  is 
greater  than  90%.  Hence,  the  total  extinction  coefficient  in  the  band 
0.92  to  I . 10  pm  is 

'total  «  l0"'  k""'- 


This  value  seems  unreasonably  small  for  total  attenuation  but  appears 
feasible  for  the  attenuation  at  the  selective  wavelength  of  1.06  pm. 


Also  note  that  Taylor  and  Yates  state  the  visual  range  (see  Figure  13)  in 
effect  at  the  time  of  their  measurements.  A  rough,  but  useful  relation¬ 
ship  that  can  be  used  to  estimate  the  total  scattering  extinction  coef¬ 
ficient,  6  =  om  +  ya ,  for  their  various  sets  of  measurements  is  [78] 


S  % 


(3.9/VM0.53/A) 


0.62V,/3 


(34) 
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where  V  is  the  visual  range  expressed  in  km  and  X  is  the  monochromatic 
beam  wavelength  in  micrometers.  The  value  of  6,  for  the  5.5  km  path 
length  is 

0,  nc  ^  4  x  ICf2  km"1 
I .  Uo 

and  for  the  16.25  km  path  length  is 

-2  -I 

0 -  n,  ^  2  x  10  km  . 

I  .(Jo 

Hence,  one  might  conclude  that  the  scattering  mechanism  was  not  as  active 
during  the  measurements  of  Taylor  and  Yates  as  it  was  for  the  measurements 
of  Curcio,  Knestrick,  Cosden,  and  Drummeter  C59]. 

If  one  removes  the  total  scattering  component  from  Taylor  and  Yate's  value 
for  Y-j-Q-j-gj  for  +he  16.25  km  path,  the  following  value  of  y '  results  where 

'I .06  =  Y total  -  B  =  4.3  x  I0'2  km'1. 

If  the  total  continuum  absorption  coefficient  is  of  the  order  of 

klj  1.5  x  I0"2  km"1 

as  discussed  previously,  then 

t|'.06  “  'total  -  B  -  kC  =  2.8  X  I0'2  km'1. 

Hence,  i+  follows  that  the  selective  molecular  absorption  coefficient  from 
Taylor  and  Yates  data  for  the  wavelength  1.06  pm  falls  in  the  range 

2.8xl02km1<k  <  6-.-3  x  10  2  km  1 

m 

when  all  attenuating  mechanisms  are  taken  into  account. 

In  section  2  of  this  report,  the  factors  that  contribute  to  the  atmos¬ 
pheric  attenuation  of  laser  energy  were  identified.  These  factors  are 
listed  in  the  table  below  along  with  numerical  values  as  given  in  section 
3.  There  is  little  doubt  that  for  all  practical  purposes  the  aerosol  is 
the  main  attenuating  mechanism  of  1.06  pm  wavelength  radiation. 
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TABLE  III  MAGNITUDES  OF  ATMOSPHERIC  EFFECTS  AT  1.06  pm 


Attenuating  Mechanism 

Attenuation 
Coefficient 
Upper  Bound 
( km“ 1 ) 

Attenuation  Loss 
Upper  Bound 
(dB/km) 

Reference 

Page  No. 
on  this 
Report 

Molecular  Scattering,  o 

a  m 

0.0008 

0.003 

[19] 

12 

Selective  Molecular  Absorp 

tion,  k 
*  m 

i-  0.065 

0.282 

1 - 1 

in 

O'* 

i _ i 

1 5 ' 

Total  Continuum  Absorp¬ 
tion,  k^ 

0.015 

0.065 

[55] 

25 

Aerosol  Extinction,  v 

a 

0.113 

0.490 

[57] 

26 

Total  Extinction, 

Y  =  o  +  k$  +  k  +  y 

m  1  m  'a 

0.2 

0.87  (clear  day,  sea  level  paths) 

Light  Rain  (1  mm/hr) 

0.2 

0.87 

[71] 

32 

Heavy  Rain  (10  mm/hr) 

1.2 

5.2 

[71] 

32 

Cloudburst  (100  mm/hr) 

6.3 

27.5 

[71] 

33 

Haze 

2.3 

10 

[63] 

27 

Fog 

58 

252 

[63] 

30 

Cloud 

132 

573 

[70] 

31 

The  results  of  this  report  indicate  that  the  total  extinction  coefficient 
at  A  =  1.06  pm  for  horizontal  paths  at  sea  level  under  good  visibility 
(>I0  km)  should  lie  in  the  range 

2  x  1 0  ^  km  *  <  v .  .  ,  <  2  x  10  *  km  * . 
tota  I 

The  attenuation  due  to  rain  over  short  paths  (<5  km)  and  at  low  rain  rates 
( < 1 0  mm/hr)  is  not  severe,  however,  cloud,  fog  and  haze  will  have  a  con¬ 
siderable  effect  on  the  transmission.  The  problem  regarding  the  trans¬ 
mission  of  1.06  'pm  vs  1.54  pm  is  still  unresolved.  The  most  definitive 
statement  that  can  be  made  at  this  time  is  that  transmission  at  these  two 
wavelengths  is  higher  than  at  most  of  the  other  laser  wavelengths  avail- 
ab  le. 

This  study  has  defined  the  following  areas  as  barrier  problems  which  need 
to  be  resolved  before  adequate  predictions  of  laser  systems  performance 
can  be  made: 

I.  Effects  of  temperature  and  pressure  on  the  half-widths  and  intensities 
of  vibrational  and  rotational  lines  are  not  adequately  known. 
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V 


TABLE  I  I i  MAGNITUDES  OF  ATMOSPHERIC  EFFECTS  AT  1.06  pm 


Attenuating  Mechanism 

Attenuation 
Coefficient 
Upper  Bound 
( km” 1 ) 

Attenuation  Loss 
Upper  Bound 
(dB/km) 

Reference 

^age  No. 
on  this 
Report 

Molecular  Scattering,  a 

3'  m 

0.0008 

0.003 

[19] 

12 

Selective  Molecular  Absorp 

i-  0.065 

0.282 

[49,51] 

15- 

tion,  k 
'  m 

Total  Continuum  Absorp- 

0.015 

0.065 

[55] 

25 

tion,  k^ 

Aerosol  Extinction,  v 

a 

0.113 

0,490 

[57] 

26 

Total  Extinction, 

0.2 

0.87  (clear  day,  sea  level  paths) 

y  =  a  +k¥+k  +  y 

m  1  m  a 

Light  Rain  ( 1  mm/hr) 

0.2 

0.87 

[71] 

32 

Heavy  Rain  (10  mm/hr) 

1.2 

5.2 

[71] 

32 

Cloudburst  (100  mm/hr) 

6.3 

27.5 

[71] 

33 

Haze 

2.3 

10 

[63] 

27 

Fog 

58 

252 

[63] 

30 

Cloud 

132 

573 

[70] 

31 

The  results  of  this  report  indicate  that  the  total  extinction  coefficient 
at  X  =  1.06  pm  for  horizontal  paths  at  sea  level  under  good  visibility 
(>I0  km)  should  lie  in  the  range 

2  x  1 0  ^  km  *  <  |  <  2  x  1 0  *  km  * . 

The  attenuation  due  to  rain  over  short  paths  ( <5  km)  and  at  low  rain  rates 
(<I0  mm/hr)  is  not  severe,  however,  cloud,  fog  and  haze  will  have  a  con¬ 
siderable  effect  on  the  transmission.  The  problem  regarding  the  trans¬ 
mission  of  1.06  'pm  vs  1.54  pm  is  still  unresolved.  The  most  definitive 
statement  that  can  be  made  at  this  time  is  that  transmission  at  these  two 
wavelengths  is  higher  than  at  most  of  the  other  laser  wavelengths  avail¬ 
able. 

This  study  has  defined  the  following  areas  as  barrier  problems  which  need 
to  be  resolved  before  adequate  predictions  of  laser  systems  performance 
can  be  made: 

I.  Effects  of  temperature  and  pressure  on  the  half-widths  and  intensities 
of  vibrational  and  rotational  lines  are  not  adequately  known. 
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2.  Vibration-rotation  interaction  in  computing  line  spectra  needs  to  be 
accounted  for  consistently. 

3.  Magnitude  of  continuum  absorption  is  unestablished  and  its  definition 
is  vague. 

4.  Contribution  of  CO2-N2I]  and  [02^2  complexes  to  the  over-all  attenua¬ 
tion  at  1.06  um  is  uncertain. 

5.  Aerosol  attenuation  measurements  versus  theoretically  calculated  aero¬ 
sol  coefficients  need  to  be  investigated  for  validity.  Controlled 
experimental  measurements  in  aerosol  absorption  work  are  seriously 
lacking,  and  research  into  the  altitude  variation  of  the  aerosol 
extinction  is  urgently  needed.  In  addition,  the  theoretical  treat¬ 
ment  of  irregular  particles  needs  further  investigation. 
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